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MULTIPLE  ELECTROMAGNETIC  SCATTERING  FROM  A  CLUSTER  OF  SPHERES.  I.  THEORY 


1.  INTRODUCTION 

Scattering  of  tight  by  molecules  is  commonly  dealt  with  through  the 

Ray leigh-Oebye  theory.  ’  This  approach  Is  known  to  be  applicable  to 

molecules  whose  effective  index  of  refraction  is  close  to  unity  and  to 

imply  a  number  of  approximations  which  may  have  rather  severe 

effects.  ^  In  this  paper,  we  face  the  problem  of  scattering  of 

electromagnetic  waves  by  adapting  to  the  method  devised  and  successfully 

used  by  Johnson^  to  calculate  the  electronic  states  of  large 

o 

molecules.  Accordingly,  we  model  a  molecule  as  a  cluster  of  spherical 
scatterers,  possibly  of  different  radii  and  complex  refractive  indexes.  A 
plane  electromagnetic  wave,  incident  to  the  cluster  undergoes  multiple 
scatterings  which  we  account  for  by  expanding  the  scattered  wave  as  a 
multi  centered  series  of  multi  poles.  The  expansion  coefficients  turn  out  to 
be  the  solutions  of  the  system  of  linear  equations,  obtained  by  expanding 
the  incident  field  in  a  series  of  multipoles  within  the  spheres  and 
imposing  the  boundary  conditions  across  the  surface.  Due  to  the  presence 
of  the  incident  plane  wave,  the  above  system  is  nonhomogeneous  so 
that  the  expansion  coefficients  are  uniquely  determined. 

The  approach  outlined  above  is  of  general  applicability  as  it  is  not 
based  on  any  assumption  -  neither  on  the  radii  and  refractive  indexes  of 
the  single  spherical  scatterers,  nor  on  the  geometry  and  size  of  the 
cluster.  The  only  approximation  required  is  the  truncation  of  the 
multipolar  expansions  in  order  to  get  a  system  of  finite  size.  The  number 
of  terms  to  be  retained  in  the  multipolar  expansions  in  order  to  get  a 
fairly  convergent  scattered  field,  as  well  as  a  number  of  related  topics, 
will  be  discussed  In  section  5. 

2  MULTIPOLAR  EXPANSIONS  OF  THE  FIELDS 

The  cluster  whose  scattering  properties  we  want  to  study  Is  composed  of 
N  nonmagnetic  spheres  whose,  centers  He  at  Ra  and  whose  radii  and  (possibly 
complex)  refractive  indexes  are  ba  and  na,  respectively.  We  refer  the 
cluster  to  a  fixed  system  of  axes  and  choose  the  direction  of  incidence  of 
the  incoming  plane  wave  through  the  direction  cosines  of  Its  wavevector. 
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A  straightforward  calculation  along  the  lines  sketched  by  Jackson3  allows 
writing  the  field  of  a  circularly  polarized  plane  wave  of  wave  vector  k  as 

1-..J  /'\1  .  , 

(la) 

(lb) 


*  £  «nU,(k)[JL(lcr)XlM(r)  *  n^lkrlXjr)] 

LM 


1B^ 

-n 


with  if  ■  +  1  according  to  the  polarization  and 

WnLM(^  =  4n1'L<?i  +  in?2)  *  (2) 

where  e?  and  e3  are  unit  vectors  orthogonal  to  k  and  to  each  other.  The 
vector  spherical  harmonics,  Xlm,  are  defined  according  to  Jackson. 9 

The  field  scattered  by  the  cluster  is  expanded  in  a  mul ticentered 
series  of  multipoles  including  only  outgoing  spherical  waves  at  infinity 


with  rQ  ■  r  -  Ra.  The  superscript  (1)  on  the  spherical  Hankel  functions  of 
the  first  kind  will  be  omitted  throughout  for  simplicity. 


As  regards  the  field  within  the  spheres,  we  remark  that  our  theory  can 
even  be  applied  to  a  cluster  of  nonhomogeneous  spheres,  provided  the  na's 
are  spherically  symmetric,  i.e.  na  »  na(ro).  Therefore,  within  the  e  -th 
sphere,  we  write 
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LM 


(4b) 


where  R*  and  S“  are  the  solutions,  regular  at  ra  »  0,  of  the  equations 
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respectively.  Of  course,  for  uniform  n^s,  R“  ■  S“  ■  jL(Kara),  with 
Kfl-  kna. 

3.  EQUATIONS  FOR  THE  COEFFICIENTS 

The  expansion  coefficients  a£lm,  B®^,  cnLM»  and  °nLM  in 
equations  (3)  and  (k)  are  uniquely  determined  by  the  boundary  conditions 
for  E  and  B  at  the  surface  of  each  of  the  spheres.  Therefore,  we  need  to 
rewrite  equations  (1)  and  (3)  In  terms  of  multipoles  centered  at  a  single 
site,  say  Ra.  This  can  be  done  by  means  of  the  appropriate  addition 
theorems  ’  which,  near  the  surface  of  the  a  -th  sphere,  l.e.  for 
ra  S  Rffl0  ■  |Rp  -  Ra|,  yield 
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and  analogous  expressions  for  1B?  and  IB^.  In  equation  ( 3-1 )  we  define 


L’M’LM 
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while  in  equation  (7),  t/“M,LM  and  z“MlLM  are  identical  to  and 

K  l'M'lm  >  respectively,  but  for  the  substitution  of  jx  to  hx  and  R^  •  0. 
The  Clebsch-Gordan  coefficients  are  defined  according  to  Rose^  and  the 
quantities 


UL'M’jLM) 
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are  the  well-known  Gaunt  integrals. ^ 

Now  we  take  the  dot  product  of  equations  (4),  (6),  and  (7)  in  turn  with 
fftY and  and  9et  the  radial  and  tangential 

components  of  the  field  at  the  surface  of  the  a-th  sphere.  Imposition  of 
the  boundary  conditions  and  integration  over  the  angles  then  yield,  for 
each  a,  l  ,  m,  six  equations,  among  which  C®LM  and  0®LM,  the 
coefficients  of  the  internal  field,  can  easily  be  eliminated.  This 
possibility  allows  getting,  for  each  a,  i,  m,  two  equations  Involving  only 
the  A's  and  B's  as  unknowns 
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where  we  define 


Rl  = 
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ct  a  a  a 


(10a) 


*£  = 


jt(kra)dT(raS“)  “  n?“dHrA(kra,) 

_ a _ a _ 

hPkrPdr~^raSP  '  naS«.dr^rah«-^kra^  r  =b 
a  a  a  o 


Pa  =  Ja  +  La 
n«wntM  ^£mLM  nl,£mLM 


Qa  -  La  +  n J* 

n  >  £m  lm  judlm  £m  lm 


(10b) 


(11a) 


(11b) 


The  system  composed  of  equations  (9a)  and  (9b)  for  all  values  of  o,  l,  m, 
completely  solves  our  scattering  problem. 

4.  THE  CROSS  SECTIONS 

Once  the  coefficients  of  the  scattered  wave,  a“lm  and  are 

known,  all  of  the  scattering  properties  of  the  cluster  can  be  easily 
calculated.  For  this  purpose,  it  is  convenient  to  express  the  scattered 
field  In  terms  of  multipoles  centered  at  a  single  point,  say  Rq,  through 
the  addition  theorem  already  used  in  the  preceding  section.  If 
rQ  -r  -  RQ,then  we  have 

^  -  EE  £  [C  lA.  (kr,)Xl1r  (r„)  ♦  L>V.  (kr0)XL.„.  (r,)l 

HIM  1  LM  L  J 
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=  5.  *,l.  <kro)?L-»’  <r«>  +  BnLW  I^L-  <r»> 


v 


and  an  analogous  expression  for  iB^*  ,  with 

a  =  EL  [a21  j°a 

+  Ba 

nCM'  hlm  l’M'lm 

a  lm  L 

HLM  L’M'  LM 

b  .  -  el  W  i°a 

+  Ba  J°a 

HL  M  nLM  U'M'LM 

01  LM  L 

r)UM  L'M'LM 
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J°a  and  L°  a  are  given  by  equations  (8a)  and  (8b)  with  *k  substituted 
for  h^.  Equation  (12)  is  valid,  provided  that  rQ%  Roa«  |Ra  -  R  1, 
i.e.,  in  the  region  outside  a  sphere  centered  at  Rq  and  including  the 
whole  cluster.  Therefore,  choosing  Rq  -  o  and  thus  letting  the  center  of 
the  expansion  (12)  coincide  with  the  center  of  symmetry  of  the  cluster,  the 
vclume  of  the  aforementioned  sphere  is  minimized.  Anyway,  the  coefficients 
an<*  »  unH*<e  those  of  the  field  scattered  by  a  single  sphere, 

depend  on  the  direction  of  the  incident  wavevector,  k.  As  a  consequence, 
all  of  the  quantities  of  interest  depend  both  on  k  and  on  the  scattered 
wavevector  k$  ■  kr,  except,  of  course,  the  scattering,  absorption,  and 
total  cross  sections,  which  depend  only  on  k.  A  straightforward 

calculation  shows,  in  fact,  that 
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1  1  nL'M’  n  l’M' 


(tot)  _  4tt2 
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(l^a) 
I2 |(l4b) 
(14c) 


Finally,  we  notice  that  the  cross  sections  depend  on  the  polarization  of 
the  incident  wave,  q  ,  as  explicitly  indicated  in  equation  (14a)  and  (14b). 

5.  DISCUSSION 

In  order  to  assist  in  discussing  both  the  physical  content  and  the  rate 
of  convergence  of  the  theory  developed  in  the  preceding  sections,  let  us 
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rewrite  the  system  of  equations  (9a)  and  (9b)  in  matrix  form: 


According  to  Waterman,1^  equation  (15)  defines  the  matrix  on  the 

left-hand  side  as  the  inverse  of  the  electromagnetic  T-matrix  for  the  whole 
cluster.  The  above  matrix  is  non-diagonal,  for  the  cluster  lacks  the  full 
spherical  symmetry;  whereas,  the  matrices  within  it  have  an  interesting 
physical  meaning  of  their  own.  The  matrices  R  and  S  are,  in  fact,  the 
direct  sum  of  the  diagonal  matrices  R  and  S  which  in  turn  form  the 
electromagnetic  T-matrix  for  the  a-th  sphere  in  the  absence  of  any  other 
scatterer.  The  presence  of  more  than  one  scatterer  in  the  cluster  is 

q  a 

accounted  for  not  only  by  the  matrices  Rp  and  S  ,  with  p  ^ a,  but  also  by 

the  matrices  H  and  K  which  couple  all  the  scatterers  to  each  other.  The 

elements  and  LM  are  shown  in  the  appendix  to  be  the  matrix 

elements,  in  the  site  and  angular  momentum  representation,  of  the  free 

space  dyadic  Green's  function.  As  the  above  quantities  appear  as  the 
coefficients  of  the  addition  theorem  we  used  in  section  *»,  this  latter, 
besides  being  a  useful  mathematical  tool,  describes  the  propagation  to  the 
site  a  of  the  spherical  vector  waves  scattered  by  the  site  0.  Therefore, 
our  previous  statement  that  multiple  scatterings  are  accounted  for  by 
expanding  the  wave  scattered  by  the  whole  cluster  in  a  mul t i centered  series 
of  multipoles  remains  fully  justified. 

The  theory  discussed  thus  far  is  based  on  general  grounds  and  requires 

no  approximation,  apart  from  the  truncation  of  the  multipolar  expansions. 

In  this  connection,  a  fairly  good  rate  of  convergence  is  expected  even  when 

the  cluster  is  not  small  in  comparison  to  the  incident  wavelength  provided 

kba  «  1  for  any  a.  Indeed  and  s“  are  known  to  decrease  rapidly  with 

increasing  l  so  that,  for  small  kba,  tf?1  and  5°  are  quite  sufficient  to 

16  17 

describe  the  scattered  wave  even  when  na  is  not  close  to  unity.  ’  ' 
Thus  the  rate  of  convergence  of  our  approach  depends  upon  the  behavior  of 
an<*  According  to  their  definitions,  the  order  of  magnitude  of 

equations  (8a)  and  (8b)  is  determined  by  the  Gaunt  integrals,  1^,  and  by 


the  spherical  Hankel  functions,  h^(kRa^).  The  1^  integrals  do  not  vanish 
when  |  si  -  L  J  <  X<  i  +  L,  but  decrease  very  rapidly  with  increasing 
X.  Thus,  although  the  imaginary  part  of  h^,  n^(kRa^j),  tends  to 

increase  when  X>  kRap,  the  eventual  effect  is  to  decrease  the  magnitude 
both  of  H and  of  with  increasing  l,  L  and  R„*.  This  behavior  was 

to  be  expected,  for  when  the  intersphere  distance  increases,  the  present 
theory  should  reduce  to  that  of  the  scattering  from  N  spherical  scatterers 
without  any  multiple  scattering  effect.  As  a  consequence,  it  is  reasonable 
to  expect  that  the  present  approach  converges  well  by  truncating  the 
multipolar  expansions  at  LM  =  3.  Since  the  order  of  the  system  (15)  is 
dM  =  2N(LM  +  I)2  -  2N,  we  should  have  -  3ON,  a  rather  high 
number  even  for  small  clusters.  However,  if  our  clusters  possess  symmetry 
properties,  as  is  the  case  for  actual  molecules,  we  can  use  group  theory  to 
get  the  system  (15)  in  factorized  form.  The  application  of  group  theory  to 
the  present  approach  to  multiple  electromagnetic  scattering  will  be  the 
subject  of  another  paper. 
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APPENDIX 


MATRIX  ELEMENTS  OF  THE  DYADIC  GREEN'S  FUNCTION 


The  free  space  propagator  for  spherical  vector  waves  (dyadic  Green's 
function)  is  the  solution  of  the  inhomogeneous  Helmholtz  equation. 

(v2  +  k2)G(r,r')  *  -4irl6(r-r')  (A-l) 

in  spherical  coordinates  and  can  be  written  with  respect  to  the  molecular 
sites  as 


fi(r.r')  * 


(A-2) 


If  we  expand  the  unit  dyadic,  1,  with  respect  to  a  spherical  basis* 

?  *  „•  £  £ (*-3) 


“jj’  J-y.'-u 


and  assume  jra  -  Rafl|  >  Neuman  expansion  of  G  is** *** 

’  4"ik?  £  hL(k|V?a6l,VLM(C.‘?aS,Jl.(kr6)Vw(fs)5-l15-u 


Now  the  addition  theorem  for  scalar  Helmholtz  harmonics  can  be  applied  to 

Vk|.r«-?j|V;»’  V ,o  «•**** 

?‘vr;>  '?££  <?„B) 

X  JL-  <kr,>\.„.  (A-M 


where  we  assumed  ra  <  Ra|j  and  thus  define 


U.<W  *  "*tk  E  'L"l"'x<L'H,‘LM)hA<kl<a6)C-»<?a6) 


*Rose,  E.M.  Multiple  Fields.  John  Wiley  &  Sons,  Inc.,  New  York,  New 
York.  1955. 

**Goertzel  G.  and  Traill,  N.  Some  Mathematical  Methods  of  Physics. 
McGraw-Hill,  New  York.  I960. 

***Nozawa,  R.  J.  Math.  Phys.  2  0966). 
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Now  we  recall  that  the  spherical  harmonics  and  the  irreducible  spherical 
tensors  are  related  through  the  equation* 


5-uWr)  SE  C(1,L,J;-u,M)t“  (r) 

so  that  equation  (A-4)  can  be  rewritten  as 

S(r  ,r')  =E  £  £  £  jL(kri)T"'pW(ri)C(l,L,J;-P,M) 
'  '6  U  JJ'  LM  L-M-  L  6  ~JL  ~B 


x6,m-  LM(Ro6)c(1»L,.‘J‘;-M.H*)jL.(kro)i;;p(ra) 
which,  through  the  position  M  -  p  »m,  M*  -  p  •m>  takes  the  final  form 

?<V:e>-£  2  E  (r„)  (A-s) 

Jm  Jm'  LL 

Equation  (A-5)  shows  that  the  quantities 

CJT',JL  =  ^  ,l^  ’  ^  ®  l*  m‘  +u ,  Lm+p  '  *  _P  ^  (A-6) 

are  just  the  matrix  elements  of  G  with  respect  to  the  irreducible  spherical 
tensors.  Moreover,  direct  comparison  of  equation  (A-6)  with  equations  (8a) 
and  (8b)  shows  that 


“l*  CM  ■  -fC  ,CC<W 


(A-7a) 


<Icm  •  ^Wcc  <«*>  •  IJWC  ,M  ™ 

Now,  since 

jc(kr>?c«<r>  -  ?«<:>  -  -Jc(kr)l"<r) 

ivxjc(kr>*LM(r)  *  Nc-<:>  ■  '[•JStJcmI".i 

an  easy  but  lengthy  calculation,  with  the  help  of  the  formulas  of  Borghese 
et  al . ,*  shows  that  H  and  X  are  the  matrix  elements  of  G  with  respect  to  M 


♦Rose,  E.M.  op  cit. 

♦♦Borghese,  F.,  Denti,  P.,  Toscano,  G.,  and  Sindoni,  0.  I.  J.  Math.  Phys. 
21,  2751*  (1980). 
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and  N.  The  other  matrix  elements  of  6  do  not  appear  in  the  present  work 
because  we  deal  with  solenoidai  fields  which  require  H  and  N  only  for  their 
description.  Finally,  we  notice  that  the  above  procedure  also  allows  the 
definition  of  1M  and  u  as  the  matrix  elements  of  G  with  respect 

to  H  and  N.  It  is,  in  fact,  sufficient  to  assume  ra  >  Rftj  and  consequently 
substitute  in  G,  j^(kR0^)  to  hx(kRa^)  and  in  M  and  N,  hL,(kra)  to  JL.  (kra). 
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